Correlated atomic motion and spin-ordering in bcc 3 He 



N. Gov 1 and E. Polturak 2 

1 Department of Physics, University of Illinois at Urbana- Champaign, 
1110 Green St., Urbana 61801, IL, USA 
2 Physics Department, Technion- Israel Institute of Technology, 
Haifa 32000, Israel 

We propose a new way to treat the ordering of nuclear magnetism of solid 3 He. We argue that the 
magnetic interaction arises indirectly as a consequence of correlated zero-point motion of the ions. 
This motion lowers the energy of the ground state, and results in a coherent state of oscillating 
electric dipoles. Distortion of the electronic wavefunctions leads to hyperfine magnetic interactions 
with the nuclear spin. Our model describes both the modification of the phonon spectra, localized 
modes ("vacancies") and the nuclear magnetic ordering of bcc 3 He using a single parameter, the 
dipolar interaction energy Eq. The model yields correctly both the u2d2 symmetry of the ordered 
phase and the volume dependence of the magnetic interaction. We calculate the magnetic excitations 
in the u2d2 phase, which compare well with the measured specific-heat, free energy and entropy. 
We further give a description of the nature of the High-Field phase and the phase-transition lines 
as a function of magnetic field and volume. 
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I. INTRODUCTION 

The spin-ordered phase of bcc 3 He presents a dif- 
ficult challenge to accurate theoretical descriptionEm. 
The main problem is to explain why the transition 
temperatureB of 10 _3 K is two orders of magnitude larger 
than the nuclear dipolar interaction ~1CP 8 K. The preva- 
lent description is in terms of atomic exchange cy- 
cles involving several atoms (Multiple Spin Exchange 
model). Successively higher order cycles produce com- 
peting ferro and antiferromagnetic interactionsuFor qual- 
itative description within this Multiple Spin Exchange 
model (MSE), one needs up to six or more atoms ex- 
changing places, which becomes difficult to calculate and 
seems a priori less probabled'EL n In fact, it is not sure 
that this description converges&ETEl for larger exchange 
cycles, i.e. may not be the spin-ordering mechanism. 
This conceptual problem, despite the relative successes 
of the MSE pictureO, therefore raises the interest in an- 
other approachEl 

We propose that magnetic ordering of 3 He is a result of 
correlations in the zero point atomic motion, which how- 
ever does not require exchange of atomalJ. IChese corre- 
lations lower the energy of the ground statetL3 and result 
in a coherent state of oscillating electric dipoles, which 
also modify the phonon spectrum. We show here that 
the zero-point motion of the atoms produces an oscillat- 
ing magnetic polarization of the electronic cloud which 
interacts with the nuclear spin in each atom. The nu- 
clear spins interact with each other indirectly, through 
the electric dipolar interaction between the atoms. In ef- 
fect, this interaction is hyper-fine like, and has the right 
order of magnitude to be related to the spin-ordering 
transition. In addition, we show that the volume depen- 
dence of the on-site magnetic interaction is in agreement 
with experimental data, and that the space symmetry of 



the proposed interaction leads naturally to the distinct 
u2d2 antiferromagnetic phase. Our scenario is somewhat 
similar to the mechanism of hyperfine enhanced nuclear 
magnetism that occurs in certain insulatorstJ. 

The plan of the paper is as follows: In section II we es- 
tablish the applicability of our model of zero-point elec- 
tric dipoles in bcc 3 He. In section III we describe the 
hyperfine magnetic interactions and their volume depen- 
dence. In section IV we show how the indirect dipolar 
interactions induced by the on-site hyperfine interactions 
lead to the u2d2 symmetry of the ordered phase. In sec- 
tion V we calculate some properties of the u2d2 phase 
and compare them to the experimental data. In section 
VI we give a description of the nature of the high-field 
phase (HFP) and calculate its properties. 

II. COHERENT DIPOLES IN THE BCC 
GROUND-STATE AND TRANSVERSE 
PHONONS 

At temperatures which are high compared to the mag- 
netic interactions (T3>lmK), the atomic motion in bcc 
3 He can be treated in thp_&ame way as in 4 HeO. 

We begin by observingMo that in contrast to the situ- 
ation prevailing in gaseous He, the interatomic potential 
is highly anisotropic in the bcc crystal. The potential well 
is both very wide and anharmonic (double- well) along the 
major axes (100,010,001). This means that the atomic 
wavefunctions will be particularly extended in these di- 
rections. In order to reduce the repulsive overlap energy, 
the many-bcidy wavefunction has to include "dynamic- 
correlations"uH The current treatment using variational 
wavefunctionliil incorporates these correlations using a 
Jastrow-type wavefunction that has a two-body or higher 
order terms to represent the short-range hard-core repul- 
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sion. Using this variational wavefunction, the phonon 
spectrum is calculated through a Self-Consistent Har- 
monic (SCH) calculations. This gives a satisfactory de- 
scription of most phonon branches, except for the trans- 
verse Ti(llO) phonon, which experimentally is found to 
be much softer than calculated. Introducing cubic terms 
into the calculation to represent the anharmonicity, soft- 
ens the calculated phonon spectra overall, which ruins 
the good agreement with the phonon branches other than 
Ti(llO). This SCH+C method also breaks the delf- 
consistency of the Harmonic procedure and predicts an 
anomalous upward slope of the lowest transverse Tuo 
mode, which is not observed. 

We would like to describe the anharmonicity in a dif- 
ferent way; we treat the local motion of the atom in its 
potential well as a "local mode". The harmonic descrip- 
tion of the crystal potential misses the low- lying vibration 
of the atoms due to the shallow (double-well) potential 
along the major axes, which is the mode responsible for 
the strong dynamic correlations in the many-body wave- 
function. We therefore treat the low-energy local excita- 
tion along the major axes, and the harmonic modes as 
two independent degrees of freedom. The collective mo- 
tion of the atoms is assumed to be well described by the 
Self Consistent Harmonic calculation. In this procedure, 
we take the SCH spectrum from previous workscJ, as the 
approximate harmonic description. The "local modes" 
can be [Shown to couple only to the Ti(llO) transverse 
phonoriiil. The coupling is described in terms of an elec- 
tric dipolar matrix element, where this zero-point dipolc 
moment arises from a correlated atomic motion along the 
major axes. 

The origin of a zero point electric dipole moment is the 
highly directional zero point oscillation of the nucleus, 
which breaks the rotational symmetry of the nuclear po- 
sition relative to the electronic cloud. Going beyond the 
Born-Oppenhcimer approximation, the electronic wave 
function aquires a p-component in addition to the s- 
orbital. The first order correction to the electronic 
ergy due to relative nuclear-electronic fluctuations 
AE 



ifl: 



(m/M n )(a/D)E s 



5K, where m is the elec- 



tron mass, M n the nuclear mass, a ~ 0.25A is the Bohr 
radius, D ~ 0.8A is the spread of the excited nuclear 
wavefunction, and E sp the excitation energy of the He 
atom from the s to the p shell. Since the crystalline po- 
tential is anisotropic, we can take this zero point motion 
to be highly directional (i.e. a completely off-diagonal 
matrix element). We assume that AE is the energy as- 
sociated with the mixing of the s and p electronic levels: 
\i>) ~ \s) + A \p), i.e. A ~ y/AE/E sp ~ 1CT 2 . This lo- 
cal polarization results .in zero-point dipolar interactions 
between the atomstilO. Taking these interactions into 
account, the ground state of the crystal can be lowered 
by having the nuclear motions correlated with respect to 
the relative phase of this zero-point oscillatory motion. 
The dipolar interaction energy Eq, of the order of AE, 
is the energy we associate with the correlated zero-point 



oscillations of the atoms along the major axes. 

In our scenario, the ground-state of the crystal in which 
the zero-point motion of the atoms is correlated may be 
described as a global state of quantum resonance be- 
tween the two degenerate configurations shown in Fig.l, 
each of which minimizes the instantaneous dipolar in- 
teraction energyEil for dipoles along the major axes (i.e. 
A = x, y, i) 



Edip — — ^ • Mi 



3 cos 2 {p. ■ (r - Ti)) - 1 



(1) 



where ro is the position of an atom and the summation 
Yi runs over the entire lattice. The atomic electric dipolc 
moment |/z| induced by the coherent zero-point motion is 
proportional to A through 



|m| =e(iP\x\ij) ~2e\(s\x\p) 



(2) 



Due to the lower symmetry of the dipolar array (Fig.l) 
compared to the crystal symmetry, a coupling between 
the harmonic phonons and the dipolar modulation ex- 
ists only along the (110) directiontil. The excitations 
of these dipoles are periodic modulations of the relative 
phases of the dipoles, with respect to the ground state ar- 
ray (Fig.l). We found that the dipolar interactions van- 
ish at the edge of the Brillouin zone, so this transverse 
phonon has the energy Eq at this point. The hybridiza- 
tion of the dipolar excitations with the phonons results 
in a factor 2 softening of the transverse Ti(110) phonon 
compared to the harmonic calculation. The experimental 
data give a softening by a factor of ~1.7, in close agree- 
ment. In addition, in our model there appears a local- 
ized (i.e. dispersionless) excitation branch of energy 2Eq, 
which is involved in mass diffusion and contributes to the 
specific- heatO. This excitation is a quantum analogue-of 
a point defect ("vacancy") and has similar featuresllil. 
The specific heat of the localized mode is consistent with 
the experimental contribution attributed to point defects 
to the total heat capacity, and in fact resolves previous 
inconsistencies related to the [properties of the point de- 
fects in this solid (vacancies)liJ. Let us emphasize that 
the additional modes are a consequence of [introducing 
an additional degree of freedom to the soliaE3, which is 
the relative phase of the internal mixing of the |s) and 
\p) states. 

In the case of bec 4 He the bare local mode energy Eq 
was measured directly by NMR of a 3 He impurity, and 
the phonon spectrum was measured using neutron scat- 
tering. Our calculation is consistent with both sets of 
datalil There are no neutron scattering measurements 
of the phonon spectrum of bec 3 He, but sound velocity 
data indicates that the slope of the Ti (1 1 0) phonon is 
about half of the SCH calculationHE3. We predict that 
this ratio should indeed be l/2p-|and we therefore take 
half the energy of the calculated^ SCH Ti(110) phonon 
at the edge of the Brillouin zone to be the bare dipole-flip 
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energy Bo(— 5K at V=21.5cm 3 /mole) in solid 3 He. Ac- 
cording to our model, the energy of the localized mode 
involved in thermally activated self-diffusion should be 
2Eq ~ 10K. This value is in excellent agreement with the 
activation energy measured by. x-ray diffraction, ultra- 
sonics and NMR experimentsEj at V=21.5cm 3 /mole. A 
similar activation energy is also obtained from the excess 
specific heatEH, and pressure measurementsEa 

We therefore establish the likely occurance of coherent 
zero-point dipoles in the ground-state of bcc 3 He. The 
volume dependence of the dipolar energy (i.e. the 'va- 
cancy' activation energy 2E(A,is. shown in Fig. 2a from 
various experimental sourcescJ'tS. Since the excitation 
energy of the atom inside the crystal potential increases 
with decreasing volume, we indeed expect the energy 2Eq 
to behave as in Fig. 2a. The mixing parameter A depends 
on pressure through y/VEo, where V is the molar volume. 
This behavior arises from the definition of Edi P 
and the factor |ro — | 3 in this equation is proportional 
to the molar volume V. We plot the result in Fig. 2b, 
where we see that A increases slowly-with decreasing mo- 
lar volume V. Future x-ray studies^ in bcc 3 He may be 
able to check our prediction for the relation between the 
activation energy (2Eq, Fig. 2a) and the Ti(llO) phonon 
as a function of the molar volume. 



III. MAGNETIC INTERACTIONS 

We now describe the magnetic interaction arising 
within our model. The lowest \p) level of the He atom has 
the electrons, in a spin S = 1 state due to strong exchange 
interactionc.3, of the order of 0.25eV. The zero-point elec- 
tric dipole interactions, arising from the relaxation of the 
Born-Oppenheimer approximation, are also accompanied 
by electronic spin interactions due to magnetic-dipole 
transitions. The size of the electric dipole matrix is given 
by the electric dipolar interaction energy Eo = 2\Edi P \ 
([j]). The size of the magnetic-dipole matrix is related to 
the electric dipole matrix byc-3 

Emagtrans 

~ E ■ {h/M a c s )/ (x) (3) 

where M a is the atomic mass, c s is the velocity of sound 
of the Tno phonon. We find that h/mc s ~ (x) ~ 0.8A 
so that both transition matrices are of the same or- 
der. We therefore have oscillating zero-point electric and 
magnetic dipoles, with an electric dipole moment of size 
~ 2eA (x) and a magnetic dipole moment of size ~ 2An e , 
where n e is the electron's magnetic moment. 

In addition, the atomic p-state, S = 1 level is split into 
3 sublevels with quantum numbers: J = L + S = 0, 1, 2, 
due to spin-orbit couplingc-3, which is of the order of 
~ 1.5K. In our particular sp-hybridized states, which are 
not pure states of the quantum number J, there is no 
spin-orbit coupling, as the orbital current is zero. Rather, 
these states are split as a consequence of the magnetic 



dipolar interaction between the two electrons. The mag- 
netic interaction will tend to align the total electronic 
spin with the symmetry axis of the excited \p) state, i.e. 
along the normal axes (Fig. 3). This follows from a simple 
calculation of the magnetic dipolar interaction of the two 
electrons, in which the |s) electron is centered on average 
in the nucleus, while the \p) electron is spread over the 
two lobes (Fig. 3). We thus estimate the aligning dipo- 
lar interaction energy at E es ~ 1.4K, of the order of the 
spin-orbit coupling. 

In the ground state, the zero-point \p) electrons will 
therefore have an oscillating magnetic moment N e ~ 
2An e , where n e is the magnetic moment of an electron. 
Because part of the electronic magnetic moment is now 
associated with the \p) state, there appears a net un- 
canceled moment of equal size in the |s) component of the 
electronic wavefunction. In 4 He the nuclear spin is zero, 
and the zero-point magnetic polarization of the electrons 
has no effect. In contrast, in 3 He the nuclear magnetic 
spin 1 = 1/2 will interact with the oscillating electronic 
magnetic moment, mainly due to the contact term of the 
|s) electron at the nucleus. The magnetic interaction is 
of the hyper-fine typecil, and the energy associated with 
it, E mag , is given by 

£ TOafl = /-yN e -N„<5(r)\ (4) 

where N n is the effective nuclear magnetic moment, re- 
duced empiricallyO by a factor of_CL83 due to zero-point 
fluctuations of the.aiuclear spinsEj. Using experimen- 
tal hyper-fine datacj, we find the maximum value of 
Emag/ks — 4.5mK (for V=24 cm 3 /mole). This energy 
is much larger than the direct nuclear dipole-dipole in- 
teraction (~ 0.1/xK), and is of the right magnitude to 
explain the high transition temperature of nuclear order- 
ing in bcc 3 He. A diagram of the energy scales involved 
in the magnetic interactions is shown in Fig. 4, where we 
show the cascade of the different interactions that ulti- 
mately align the electronic and nuclear spins, with the 
zero-point excitations of the p-state. 

Before proceeding to describe the nuclear spin ordering 
phenomenon, and the effect of the magnetic interactions 
on the zero-point electric dipoles, we present a quantita- 
tive calculation of the volume dependence of the magnetic 
interaction at T=0. It is known empirically that these 
interactions have a strong volume dependenceEl, usually 
described as: E mag oc V 18 . In our model, the magnetic 
energy (||) will change with pressure due to changes of the 
electronic magnetic polarization N e . First, the pressure 
changes the admixture of the \p) orbital in the ground 
state wavefunction through changes of the mixing 
parameter A. Second, as the solid is compressed, the 3 
sub-levels (J = 0, 1, 2) broaden into partially overlapping 
bands. This overlap leads to further reduction of the net 
alignment of the electronic magnetic moment N e with 
the orbital polarization (Fig. 3). 

The reduction in N e due to the broadening of the sub- 
levels with pressure can be approximated as linearly pro- 
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portional to the overlap integral of the three spin-orbit 
sub-levels 



(N e (V))~(l-F(V))\(V)N e 



(5) 



where the overlap integral F(V) = J (j>* 1 (E)cj)„ 1 (E)dE, 
and <j)*(E),(f>—i(E) are taken as normalized Gaussians 
centered at the energies of the spin-orbit levels, which 
are ~ 1.5K apartnJ. The effective energy width 7 of 
these Gaussians (c/>*(E) and (j>_i(E)) is estimated using 
a simple band calculation, as the change of the Coulomb 
energy of the \p) electrons due to .overlapping \p) wave- 
functions on neighboring He nuclei 2 ^! 



1{V) = x 2 J2 



r ^(r)^i(r)d 3 r 



(6) 



where ipQ (r) , ipi (r) are the p-state wavefunctions of the 
electrons around the central atom and its i neighbors. 
In this calculation we also take into account the large 
spatial spraad of the atomic position inside the wide po- 
tential well", of approximately ±0.8A. We find that 7 
ranges from ^OK for V=24 cm 3 /mole to 5.5K for V=19 
cm 3 /mole. The overlap factor 1 — F(V) is therefore 
very sensitive to volume, changing from 1 at V~ 24.3 
cm 3 /mole to ~ 0.01 at V=19 cm 3 /mole. One can see 
that as the volume decreases the broadening of the bands 
increases, thereby decreasing the net magnetic polariza- 
tion of the electronic cloud. 

The strength of the magnetic interaction should be 
proportional to the measured Curie- Weiss temperature 
9. In Fig. 5 we compare the normalized magnetic inteiu 



action energy E„ 



with the normalized measure 



9. We find that the volume dependence of E mag agrees 
very well with that of 9. Since in our model there is only 
one magnetic interaction, there is a single energy scale 
that determines the volume dependence of the transition 
temperature and the e2 coefficient of the specific-heatEl 
All these quantities have the same volume dependence, 
a simpler situation than that ijawflving a competition 
among many exchange processes™! 



IV. COMPOSITE QUANTUM RESONANCE OF 
THE ELECTRONIC AND MAGNETIC 
SUBSYSTEMS 

We now consider the effects of the magnetic interac- 
tions on the coherent-state of the oscillating electronic 
dipoles. We begin with the high-temperature param- 
agnetic phase (PP), where the nuclear spins are ran- 
domly oriented. The existence of the hyper-fine splitting 
Emag @) means that the quantum resonance state of the 
atoms, oscillating between the two configurations of Fig. 1 
is broken. The resonance frequency of each site will be 
slightly different from each other, randomly shifted by 
— E mag /h, along each of the 3 major axes. The result 
will be that the zero-point electric dipoles will slowly drift 
out of resonance, breaking the long-range order of Fig.l. 



It is possible to restore the degeneracy of the overall 
ground state, and hence the quantum resonance condi- 
tion of the zero-point dipolar oscillations. This can be 
done if the nuclear spins become ordered as well, in a 
spatial configuration which has the same overall magnetic 
interaction energy with each of the two states of Fig. 1. 
We first note that the total on-site magnetic interaction 
is isotropic with respect to the direction of the nuclear 
spin. This is due to the simultaneous interaction of the 
nuclear spin with the zero-point electronic spin which has 
equal components oscillating along all 3 major axes. 

The electric dipolar interactions in the ground-state 
of Fig.l are such that each simple cubic sublattice of 
the bcc has no net interaction with the other sublattice. 
We therefore look for the possible static arrangements of 
the nuclear spins that fulfill the resonance requirement 
on each sublattice independently. Such arrangements 
are shown in Fig. 6. In these arrangements there is an 
equal number of atoms with electronic and nuclear spins 
aligned (and anti-aligned) in both degenerate configura- 
tions of the electric dipoles (Fig.l). This arrangement 
also preserves the overall time-reversal symmetry of the 
system at zero field, i.e. no net magnetic polarization. 
Since these configurations exist on both sublattices, we 
end up with the u2d2 arrangements which is the symme- 
try of the ordered nuclear phaseESEl The only constraint 
so far due to the quantum resonance is the u2d2-type 
symmetry (i.e. also u3d3 etc.). We would like now to 
describe in more detail the nature of the coherent state 
combining the zero-point electronic dipoles with the or- 
dered nuclear spins. 

The random on-site magnetic interaction can be 
thought of as a fliping of the dipoles with respect to the 
ground-state configuration of Fig.l. The average density, 
or probability, of a dipolc-flip is given by the ratio of the 
magnetic energy E mag to the electric dipole-flip excita- 
tion energy Eq. The ground-state occupation of these 
flipped dipoles due to the magnetic interactions, means 
that there is a change in the zero-point occupation of the 
p electrons A. The p-state mixing which is associated 
with these m agnetically -induced flipped dipoles, is given 
by: A m ~ \y/E mag /Eo ~ 10~ 2 A. This small polarization 
corresponds to an electric dipole moment: jj, m ~ 10~ 2 /i, 
where fi is the full zero-point electric dipole moment (j^). 
Since the magnetically-induced electric dipole-moments 
fj, m are interacting with each other, we are looking for 
configurations (of the type shown in Fig. 6) that mini- 
mize this dipolar interaction energy (Q) along the major 
axes (i.e. ft = x, y, z) 



E 



magdipole 



Mm,0 ' Mm,i 

i#0 



3 COS 2 (/2 m • (r - !•;)) - 1 



(7) 



The configuration of Fig. 6a has the lowest dipolar en- 
ergy given by (|7]), and will therefore be chosen as the 
ground-state of the ordered phase. In this configuration, 
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the nuclear spins are aligned along a vector d which is 
orthogonal to the u2d2 modulation axis 1 (here along the 
(100) axis). One can see that the total dipolar interaction 
energy will be invariant under any rotation of the vec- 
tor d around 1. We therefore naturally obtain the state 
where the vector 1 is along one of the major axes, as is 
found experimentally from the domain structureO. We 
point-out that the u2d2 symmetry of the ordered phase 
results from symmetry and energy considerations which 
are independent of any quantitative parameters. This 
is in contrast to the MSE theory where the symmetry, ( 
determined by competition between many parameters!!! 

In our picture, the interaction between the nuclear 
spins is indirect, being mediated by the electronic sub- 
system. This picture is somewhat similar to the hy- 
pcrfinc enhanced nuclear magnetism observed in certain 
insulatorsEj. There, the nuclear polarization is effectively 
enhanced by an electronic spin that is partially aligned by 
a strong hyperfine interaction of a charged rare-earth ion. 
In our case, it is a quantum resonance of the electronic 
system that results in a coherent phase of oscillating elec- 
tric and magnetic moments. 

Experimentally, there is a marked difference between 
the hep and bcc phases, in the sense that there is spin- 
ordering in the bcc phase, as .apposed to the lack of such 
a transition in the hep phased. This result is a natural 
consequence of our model, since we expect the triangular 
symmetry to frustrate any long-range dipolar orderingtil 
of the type shown in Fig.l. Indeed there is no evidence 
for long-rangej-aero-point electric dipole correlations in 
the hep phaseLL2l, in the form of large discrepancies be- 
tween the SCH calculations of the phonon spectrum and 
the measured phonons. We therefore do not expect a 
priori any -Long-range ordering of the nuclear spins in the 
hep phaset3. In contrast, within the conventional MSE 
model, the absence of ordering of the spin system well be- 
low its Curie- Weiss temperature is not easily explained. 



V. MAGNETIC EXCITATIONS AND ENERGY 
OF THE U2D2 PHASE 

Due to the large ratio between Eq and E mag , the slow 
zero-point oscillations of the magnetically-excited elec- 
tric dipoles /i m can be decoupled from the much faster 
zero-point oscillation of the electric dipoles with natural 
frequency Eq/%, The new coherent state can therefore 
be described as a tensor product of the coherent states 
of the two decoupled components of the oscillating elec- 
tric dipoles. A fast component with natural frequency 
Eq/H described by the Bose creation/annihilation op- 
erators ak\a,k representing approximately the full elec- 
tric dipole moment (||), and a slow component described 
by the Bose creation/annihilation operators ,b k rep- 
resenting the electric dipole moment [i m induced by the 
on-site hyperfine magnetic interaction. This decoupling 
of the magnetic and phononic quasi-particles is quite 



naturalB, due to the difference of 3 orders of magnitude 
in the energy. 

If we treat the magnetically- induced interactions as de- 
coupled from the fast oscillating component, the total 
effective Hamiltonian for the new coherent state is 

H = Hp] lon + H mag (8) 

where H p h on is given in Eq.5 of Ref. [0], and similarly 

Hrnag — ^ ^ {^-^jnagdipole X mag (k) ) b^bfc 
k 

+ Y,Xmag(k)(bibl k +b k b_ k ) (9) 



where b\,,bk are Bose creation/annihilation operators 
that represent a local flip of /i m out of the ground- 
state (u2d2) arrangement, which behaves as a local-mode 
with bare energy @j E magdlp = 2E magdlpole ~ E mag , 
i.e. X mag (0) — —EjnapMp/2. The dipolar interaction 



X mag (k) is given by <M 



Xmag (k) — ^ Cm,0 ' fm,i 

x exp [27rik • (r — rj] 



3 cos 2 (ft m ■ (r - r,)) - 1 



r - r L 



(10) 



with k being the wavevector characterizing the modula- 
tion of the array of the dipoles fi m and nuclear spins, 
arranged in u2d2 configuration Fig. 6a. Since the relative 
phase between the magnetically-induced electric dipole 
moments /i m is determined by the relative alignment of 
the nuclear spins, there is a one-to-one correspondence 
between the modulation of the nuclear spins relative to 
the u2d2 arrangement, and modulation of the relative 
phase of the electric dipoles (i m . The local-mode of en- 
ergy E mag di P therefore corresponds also to a nuclear spin- 
flip out of the u2d2 ground-state arrangement. 

The Hamiltonian H mag (||) which describes the effec- 
tive interaction between the magnetically induced dipoles 

be diagonalized using the Bogoliubov transformationE3 
(3 k = u(k)b k + i>(fc)fo^_fc. The two functions u(k) and v(k) 
,2(U\ _ l ( E magdip +x mas (k) ! ^\ n< 2, 
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+ i \y{k) 



— 1 ] , where the modulations of the 



are given by: u (k) 

1 / E magdlp +X mag (k) 

2 ^ E sm (k) 

magnetic Hamiltonian H mag describe the energy spec- 
trum of the spin- waves. 

^flir (^0 — \J Ematjdip {Ernagdip ~i~ 2X mag {k^) (H) 

The result of solving pthe. effective Hamiltonian (^) is a 
coherent ground-stateaO 



l*o> 



J|cxp 



\vac) 



(12) 
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A modulation of the relative phases 
of the magnetically-induced electric dipoles fj, m with a 
wavevector k, represents a spin- wave since the nuclear 
spins have the same relative phase. We see that the elec- 
tric dipolar interaction introduces an effective restoring 
force between the nuclear spins, that results in the above 
dispersion relation. 

In Fig. 7 we plot the dispersion relations of the spin 
waves obtained from Eq. (ll|), along various directions 
of the crystal. We find that the velocities of these waves 
at long wavelengths are in the range 6.0 — lO.Ocm/sec 
(for V~ 24cm 3 /mole). These values are similaii-to those 
extracted from melting pressure measurenieiitsLJ, where 
a value of 8.4 ± 0.4cm/sec was found. There is a soft ex- 
citation in the (100) direction (c sw ~ 3cm/sec), and also 
an unstable excitation with an imaginary dispersion re- 
lation in the (010) direction, along which we therefore do 
not expect a propagating spin- wave. Along this direction 
the modulation in the nuclear spins breaks the quantum- 
resonance condition, which states that the total magnetic 
interaction energy is the same for both configurations of 
the dipoles of Fig.l. This modulation is therefore pro- 
hibited within the quantum resonance restriction. We 
note that this calculation was for the nuclear spins fully 
aligned with the z-axis. I— . 

Mean-field calculations of the MSE modelEj have 
also resulted in anisotropic spin- waves, though there 
is at present no experimental data on the directional 
anisotropy of the spin-waves to compare with theory. 

In addition to Bosonic spin-waves there can also be 
local dipole-flips of the magnetically-induced electric 
dipoles fx m . These 7r-shifted, i.e antisymmetric, dipole 



excitations are treated as Fcrmions, just a: 
viously for the full electric dipole momenir 2 
spectrum of these excitations is given b^ 1 



Ef m {k) = yj 



E 2 

magdip 



E sw {k) 2 



escribed pre- 
. The energy 



(13) 



The dispersion relation of these optic-like excitations is 
also plotted in Fig. 7. 

We can now calculate the contribution of all these 
excitations to the specific-heat. In the calculation, we 
included the individual contributions of each mode (in 
a spherical phase-space volume), and then simply aver- 
aged. Each type of excitation was treated with using the 
appropriate (Bose or Fermi) statistics, and the details 
are described in the Appendix ( |l§| ) . The agreement with 
the (limited) experimental data for the volume V=24.13 
cm 3 /moleEl is fairly good (Fig. 8). At this molar volume, 
we have E mag di P — 2.6mK. We note that the contribution 
of the optic-like modes is important at higher tempera- 
tures, and that the agreement between the data and the 
simple T 3 law is somewhat fortuitous. .—Similar conclu- 
sions were reached using the MSE modeO. In particular, 
our calculation has the right slope as the experimental 
data. 

For the specific volume V ~ 24cm 3 J4pole, experimen- 
tal data on the ground-state energyE3E3 gives E mag /2 ~ 



1.1 ± O.lmK, while our calculation (^,^|) gives E mag /2 ~ 
1.3mK. Plotting the free energies for the two phases 
(fl8|),(|il|) in Fig. 9, we find that they give a transition tent, 
perature: Tn ~ 1.25mK, close to the measured resultEZl 
of ~ 1.05mK. In this calculation, we took the empirical 
value of E mag /2 = l.lmK, for the constant aligning en- 
ergy appearing in (p"8|) , (|l6|) . In addition, the experimen- 
tal data in Fig. 9 was shifted up by O.lmK, which is within 
the experimental uncertainty. We find that the flat free 
energy below the transition temperature fits the free en- 
ergy of the ordered phase ( |l8| ) . It should be pointed out 
that this calculation, which is based on independent spin- 
waves, is valid only at low temperatures, compared to 
E m ag- We see that the energy difference between differ- 
ent spin arrangements which fit the symmetry constraint 
(such as u3d3 etca), is only through the relatively small 
term AEdi P in the free-energy (|l7|). 

In Fig. 10 we show the entropies of the two phases as 
a function of temperature, and indicate the jump at the 
approximate transition temperature. We obtain a latent 
heat of AST C ~ 100/iK per atom, while experimentally 
it is found to beO ~ 250/iK per atom. This degree of 
agreement is quite reasonable considering the crudeness 
of the spin-wave approximation at temperatures where 
the entropy is already ~ 0.35 ln(2), and the independent 
spin- wave picture is questionable. 

Finally we describe the magnetic response of the u2d2 
phase. The zero-temperature susceptibility of an anti- 
ferromagiifit with the field orthogonal to the magnetic 



momentsEI, is given by: x(0) = y? n /E mag , where E mag /2 
is the aligning energy, and we use the effective nuclear 
moment, reduced by a factor of 0.83, due to the quantum 
fluctuationsO. The calculated susceptibility is smaller by 
~ 15% than the measured valueEj, and is therefore in 
reasonable agreement. 

The zero-field antiferromagnetic resonance frequencj 
measured by NMRo, is given by: ti 2 ^ = AE d E mag , 
where Ed is the direct nuclear magnetic dipolar interac- 
tion. To calculate the temperature dependence of this 
quantity, we calculate (Appendix) the reduction in the 
ordered magnetic moment due to thermal occupation of 
excited statesEJ. Since Ed depends on the ordered nuclear 
magnetic moment N n squared, we have that Qo oc N n . In 
Fig. 11 we plot the normalized temperature dependence 
of Qq compared to the thermally reduced magnetic mo- 
ment: (1 - AM) (Eq.Q). We find that the reduction 
in the frequency is well described by our calculation, up 
to ~ 0.7mK, above which the experimental data have a 
steeper decline. We further note that NMR experiments 
measure the resonance of the stationary nuclear spin, and 
do not reflect the fast-oscillating [u> — Eq/K) zero-point 
electronic magnetic moments N e . 

Our final comment about the u2d2_Djhase is concern- 
ing its domain structure. It was founcOuj that a crystal 
grown in a narrow channel will have one large magnetic 
domain which will be oriented such that the magnetic 
planes are parallel with the channel axis and surface. 
This observation can be explained by our model in the 



G 



following way. On the channel surface the coherent mo- 
tion of the atoms is inhibited in the direction normal to 
the wall. The only remaining coherent motion is parallel 
to the wall surface, and so will therefore be the magnetic 
ordering. At the nucleation site, when the crystal grows 
into the free liquid, it will have the three domains along 
the three orthogonal major axes. The relative sizes of 
the domains will be-determined by random roughness at 
the nucleation sitecil. 



VI. HIGH-FIELD PHASE 

Above a critical external magnetic field H c \ there is 
a first-order transition tOr-a, high-field phase (HFP). The 
NMR data for this phasec2l indicates a phase with cubic 
symmetry and large magnetization: N z > 0.6 (in units 
of the saturation magnetizations) . The nuclear spin 
arrangement infered for this phase is shown in Fig. 12. 
Since there is now global nuclear magnetization along 
the field direction (z-axis in Fig. 13), there is no time 
reversal symmetry and no quantum resonancea, for the 
oscillating electric dipole moments /i m along this axis. 
Along the orthogonal a;y-directions there is still quantum 
resonance for the electronic dipoles /i m 's, and hyperfine 
interactions of energy E mag . In a semi-classical picture, 
the nuclear spins are vectors with the magnetization in 
the z-direction given by their z-component, precessing 
around the z-axis with the Larmor frequency. In the 
paramagnetic-phase (PP) the xy-components of the spins 
are random, while in the HFP there is long-range order 
for these components (Fig. 13), due to the aligning effect 
of the hyperfine interactions in the xy-plane. The order- 
parameter is therefore the global phase of the coherent 
xy-components /zjjf, i.e. a complex order-parameter, and 
we expect a second-order phase transition of the type 
found for example in superfluid He, which is also in the 
3D-X.Y universality classed. Further, the specific-heat 
dataQ of the HFP has a distinct cusp at the transition 
temperature, typical of a A-type transition. 

Let us examine the field dependence of the hyperfine 
magnetic interactions in the ccy-plane, by first looking 
at the electronic polarization. With increasing external 
magnetic field the electronic spins N e in the xy-plane 
will become oriented with the external field along the z- 
direction, reducing the hyperfine interactions with the 
nuclear spins in the xy-directions (^). The electronic 
spins are oriented with the oscillating p-states of the 
atom in the xy plane, with effective energy (E es ) ~ 
1.4(1 — F)K, where the factor 1 — F takes into account 
the overlap of the different electronic spin orientations 
(||, Fig. 3), reducing the effective energy gap. The single- 
atom magnetic Hamiltonian for the electronic spins in 
the xy-plane, in an external field H is therefore 

Helec * —EesNf - HN* (14) 



with N£,N* the classical components of N e . The first 
term represents the aligning dipolar interaction between 
the s and p electrons of the atom. Solving (]lj) by 
mean- field, we get a linear magnetization: N* = H/E es , 
and full magnetization of the electronic spins N e at a 
field: H c2 ~ 26T (using (1 - F) ~ 0.6 for V- 24cm 3 , 
so that E^J2 ~ 0.4K), in rough agreement with the 
estimateoola experimental value of ~ 22T. At fields 
larger than there are no xy hyperfine interactions 
so the nuclear spins are free to fully align with the ex- 
ternal field, and there is a zero-temperature transition to 
the PP. n 
The volume dependence of H C 2 is experimentally^ 
found to have a volume dependence which is steeper than 
the critical Neel temperature Tjv of the u2d2 phase. This 
is in agreement with our formulation ( |l4| ) , as H C 2 is given 
by the steep dependence of the factor (1 — F) (Fig. 5). 
Since Tjv is roughly proportional to E mag (Fig. 9), we 

S' ct: tflSSl - (^Eqs.|,|and A(V) from 

Fig. 2a). Comparing the data£3E§c3 at V2=24 cm 3 /mole 
and Vi=22.7 cm 3 /mole we find: - 0.25, ~ 

1.2 and ^ N ^l ~ 0.3, in excellent agreement with the 

i N ( V2 ) ° 

above relation. 

If the electronic spins \x m in the xy-plane would have 
an arrangement where the dipolar interaction energy is 
negative (such as the u2d2 phase (0)), they would have 
zero-point oscillations with the resulting dipolar inter- 
action energy E mag dip (|§,|l(J)- These oscillations would 
get them out of phase with the nuclear spins, which are 
precessing with the Larmor frequency. We therfore can 
only have a cubic arrangement with zero ground-state 
dipolar interactions. In Fig. 13 we plot the dipolar in- 
teraction energy Xh f p (|Io| ) of the icy-components ^ as 
a function of modulation wavevector k along the (001)- 
axis, with respect to the ground-state array of Fig. 13. We 
find that the dipolar interaction energy is indeed zero in 
the ground-state, with a quadratic dispersion similar to 
excitations in a ferromagnet. 

We shall now describe the nuclear spins in the HFP. In 
the HFP the nuclear spins will have a stiffness against 
small perturbation from the ordered configuration of 
Fig. 13 (see Fig. 13). Perturbations of the nuclear spins 
in the ccy-plane will cost an energy X^fp (Fig. 13), while 
the z-component perturbation costs the Zeeman energy 
Ezeeman = N n H . A simple description of this state is 
therefore by the following model magnetic Hamiltonian 
of independent nuclei 

H hfp ~ -Ha z N* - E x a x N% (15) 

where H is the external field in units of energy, E x = 
E mag I ^ is the effective aligning energy of the xy- 
components compared with the random arrangement of 
the nuclear spins (see Eq.[l6]) , N* , N% are the compo- 
nents of the nuclear spin. The size of the hyperfine 
xy-interactions of the nuclear spins with the electronic 
spins will be reduced with increasing magnetic field: 



7 



ill- 



E x oc Nf ~ y/1 - H/H c2 

The energy of the nuclear spin is therefore given by 
Ehfp = -\/H 2 + E%. We plot in Fig. 14 the energy E hfp 
compared with the energy of the u2d2 phase at zero tem- 
perature P): E u2d2 = -{E mag /2) - AE dlp - x(0)# 2 



(using at V=24cm 3 E mag /2 ~ l.lmk, see Sec. V). We 
see that the HFP has lower energy than the u2d2 phase 
at fieMs-i/ > 0.7T, which is quite close to the measured 
valueEj'cl of H c i ~ 0.45T. In Fig. 14 we also plot the mag- 
netization of the u2d2 phase: N u2 m^ x(fyH, compared 
with the measured magnetizatiorHEl, and as calculated 
for the HFP ©: N hfp = -2H (H + E hfp ) /E 2 X . We 
find that the calculated magnetization is ~ 0.55 at the 
calculated critical field H c \, which is close to the mea- 
sured valucEJtl of ~ 0.6. 

We further checked the dependence of this critical field 
on the magnetic energy E mag , and we find that it has a 
weaker than lineat-dependence, in accordance with the 
experimental datacJ. Taking E mag oc Tjv an( i using the 
measuredEa values at Vi = 24.21cm 3 and V 2 — 22.45cm 3 , 
we calculate a ratio of critical fields of H cl (V\ ) / H cl {V 2 ) ~ 
3.4, compared to the measured ratio of ~ 3 (while the 
transition temperatures vary by a factor of ~ 4.2). 

We note that our description of the HFP is at best 
qualitative, since at higher magnetic fields the calculated 
magnetization Nhfp is much larger than the measured 
magnetization (Fig. 15). Using the measured magnetiza- 
tion we shall now use ([l5]) to describe the phase diagram 
of the HFP. The overall stiffness of the spins should be 
proportional to the transition temperature of HFP to PP, 
since the stiffness represents the energy cost to random- 
ize the spins and create a PP. From ( |15| ) and Fig. 13 this 
energy should be given by the geometric-mean of the 
energies appearing m 

©: E mean = y/(HN z )(E x N x ), 
where N z is the measured magnetization (Fig. 15) and 
N x = yl — N 2 . We find that a constant factor of ~ 1.3 
relates this last expression (E mean ) with T c , and is in 
excellent agreement with the measured T c (Fig. 16). 



VII. CONCLUSION 



aerogel) and in lower dimensions such as adsorbed two- 
dimensional films. 

Acknowledgements N.G. thanks Tony Leggett and 
Gordon Baym for many useful discussions and sugges- 
tions. This work was supported by the Israel Science 
Foundation, the Technion VPR fund for the Promotion of 
Research, the Fulbright Foreign Scholarship grant, NSF 
grant no. DMR-99-86199 and NSF grant no. PHY-98- 
00978. 



VIII. APPENDIX 

Here, we present details of the calculation of the free- 
energy and magnetization of the u2d2 phase. 

The ground-state energy of the nuclear spins has two 
contributions. The first is the reduction of energy com- 
pared to the random case, where the electronic dipoles 
fi m of ever age nergy E mag /2 are excited. This energy 
is equivalent to an effective magnetic field that tends to 
align the nuclear spins, and its effects already appear 
significant in the P£L The free energy of the PP, can be 
therefore written asc^l 

F P ara ^ -k B Tln [2 cosh(E mag /2k B T)} (16) 

In the ordered, coherent u2d2 phase, there is an addi- 
tional reduction in the ground-state energy of the coher- 
ent magnetically-induced el ec t rk ; dipoles [i m due to their 
negative dipolar interactiorStJ 

a \ * E sw (k^} (E mag di p -(- X rnag {k^ 

l\&div = 2^ 2 

k 



-50fj,K/atom 



(17) 



where we again calculated the contribution from each 
Bosonic mode and then averaged. In the u2d2 phase 
there are therefore coherent zero-point electronic dipoles 
Hm, which reduce the ground-state energy (p^), oscillat- 
ing with a resonance energy E mag( n p (^Jlt]) ■ The free en- 
ergy of the ordered phase can therefore be approximated 



To conclude, our model enables us to describe both 
the phonon spectra and the nuclear magnetic ordering of 
bcc 3 He using a single parameter, which is the thermal 
activation energy Eq. The model offers a consistent pic- 
ture of the nature of the magnetic phases of bcc 3 He, and 
their phase diagram. We give a quantitative description 
of the u2d2 phase, which compares well with the mea- 
sured data. Despite the more qualitative nature of our 
description of the HFP, we obtain a good description of 
the phase diagram parameters: H c \ 1 H c2 ,T c and their 
field and volume dependence. 

We are currently exploring the consequences of this 
model regarding the description of the magnetic proper- 
ties of solid 3 He in confined geometries (nano-clusters, 



F u2 d2 — ^Edi P — (E mag /2) 



+k B T 



^m(l-e X p(-£ s ™(fc)/fc B T)) 



- ^ln (1 + exp(-E fm (k)/k B T)) 



(18) 



The specific-heat C v = —Td 2 F/dT 2 and entropy S = 
dF/dT are then calculated from these free energies. 

We now calculate the reduction in the ordered mag- 
netic jjpioment due to thermal occupation of excited 
statcsEl. For the Bosonic and Fermionic excitations 
( |j"l| , p"3"| ) we have a reduction in magnetization which is 
proportional to the occupation of excited states 
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AM = J2Hn B +vl(l-n B )] 

k 

Yl [ W l n F + Vki 1 + n F)] 



(19) 



where ns, np are the Bose/Fermi occupation factors re- 
spectively, and (ufc|Hijfc), (wk,Uk) are the Bose and Fermi 
Boguliubov factorslI3. The factors (ufc,«fc) where givi 
before (section V) , while the Fermi factors are given b; 



Emagdip \ 2 ^ /, -^magdip 

E fm (k))' Vk= 2 \ ~ E fm {k) 



(20) 
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FIG. 1. The two degenerate 'antiferroelectric' dipole ar- 
rangement in the ground-state of the bcc phase. The arrows 
show the instantaneous direction of the dipoles. 
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E(spin-orbit)~1.5K 
FIG. 4. The energy tree-diagram, showing the effects on 
the electronic and nuclear subsystems due to the coherent 
"zero-point dipoles. The first effect is the excitation of A p-state 
with 5 = 1, with energy Eo- The second effect is the align- 
ment of the electronic spins with the electric dipole (E e3 ). The 
last effect is the hyperfine interaction that aligns the nuclear 
spin with the electronic spin with energy E mag (Eq.4) . 
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FIG. 2. (a) The dipelar activation energy 2Eq ifrpm exper- 
imental data: CirclesE3, squaresEil and diamondsEa. (b) The 
mixing parameter A (Eq.2) from the experimental data of (a). 



S electronic orbital 



P electronic orbital 



1.5 



=i 
05 

s— 

O) 
Q. 

E 



2 1.0 - 



I 

o 

<D 



0.5 - 





FIG. 3. The spin-orbit electronic magnetic-dipole interac- 
tion E Ba which aligns the s-shell spin in the center (dashed 
arrow) with the p-shell spin in the lobes (solid arrows). 
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FIG. 5. The normalized magnetic splitting 2E mag (Eq.4) 
as a function of volume, calculated using different experi- 
mental data sets for the activation energy Eo: dashed lineE- 1, 
dash-dotEH and solid lineE3. The solid circles i»re the normal- 
ized experimental Curie- Weiss temperature wU. 
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FIG. 6. Example of two simple u2d2 static nuclear spin 
(arrows) arrangements that maintain the quantum resonance 
of both the electric and magnetic dipoles (Fig.l). Only the 
spins on one simple-cubic sublattice of the bcc are shown, 
with the other sublattice having identical arrangement. 
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FIG. 9. The calculated free energies of the coherent u2d2 
magnetic dipoles (solid line, Eq.18), of the PP (dashed line, 
Eq.l|6^, of free spins (thin dashed line) and the experimental 
datauJ (solid circles). The inset shows the transition region. 
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FIG. 7. The spin-wave spectrum E aw (k) (Eq.ll) and 
the optic-like branches of the fermionic excitations Ef m (k) 
(Eq.13). The nuclear spins are arranged as shown in Fig. 6a, 
and the excitations are long the directions: (lOO)-solid, 
(lOl)-long dashed, (001) dash-dot, (Oll)-dashed, (lll)-dot- 
ted. 
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FIG. 10. The calculated entropies of the ordered phase 
(solid line) and of the PP (dashed-line). At the approxi- 
mate transition temperature of lmK we show the jump in 
entropy which gives a latent heat at the first-order transition 
of ~ 100/iK per atom. 
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FIG. 8. The calculated total spe- 

cific-heat (from Eq.18) (solid line), compared to experimental 
resultso at V = 24. 13cm 3 /mole (solid circles). We also plot 
the contribution of the fermionic modes (heavy dashed-line), 
and the linear spin- wave approximation (thin dashed line) . In 
this calculation we used E mag = 2.6mK. 
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FIG. 11. The temperature dependence of the calculated 
ordered nuclsar magnetic moment (Eq.19) compared with 
the measuredEj normalized antiferromagnetic resonance fre- 
quency Qq- 




FIG. 12. The High-Field Phase (HFP) arrangement of the 
nuclear spins. A ferromagnetic arrangement in each sublat- 
tice of the zy-component (dashed arrows), with zero total xy 
magnetization. The external field B is in the ^-direction. 
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FIG. 13. (a) Schematic description of the small spin per- 
turbations of the xy and z components, (b) The dipolar inter- 
action matrix Xhf P (k) between the xy electronic dipoles \i m , 
as a function of the wavevector in the (OOl)-direction (solid 
line), (lOl)-direction (dashed line) and typical spin-waves of 
the u2d2 phase (dash-dot lines). 
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FIG. 

Comparison between the calculated zero-temperature ener- 
gies of the u2d2 phase (E U 2d2 Eq.[jl| solid line), HFP (Ehf P 
Eq.^i], dashed line) as a function of magnetic field. The calcu- 
lated transition at H c \ ~ 0.7T. (b) Comparison between the 
calculated zero-temperature magnetization of the u2d2 phase 
(N U 2d2, dash-dot-Jine) , HFP (Nhf P , dashed line) and the ex- 
perimental datauJ (solid circles with guiding solid line). The 
measured transition is at H c i(exp) ~ 0.5T. 
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FIG. 16. The geometric-mean energy E mea n X 1.3 (solid 
line) compared withi-the normalized-4ransition temperature 
T/T c (0) (solid circlesu, solid squaresEj) as a function of nor- 
malized magnetic field H/H C 2- 
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FIG. 15. Compari- 

son between the calculated zero-temperature magnetization 
of the HFP (Nhfp, dashed line) and the experimental dataO 
(solid circles with solid line). 
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